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System and Method for Multi-access MIMO Channels with 
Feedback Capacity Constraint 

FIELD OF THE INVENTION 

[0001] The present invention relates to wireless communication systems. 
More specifically, the invention relates to optimizing the resource allocations of 
multiple users in a Multiple-Input, Multiple-Output (MIMO) system under a 
feedback capacity constraint. 

BACKGROUND OF THE INVENTION 

[0002] The growing popularity of wireless communication for voice, data, and 
video applications has created a need for larger wireless capacity. 

[0003] Current wireless communication systems are based on the Single- 
Input, Single-Output (SISO) channel where a single transmitter communicates 
with a single receiver over the channel. In one common arrangement, many 
mobile devices, such as cell phones, may be in communication with many base 
stations that are connected to a wired communication system. Either the mobile 
device or the base station selects a channel over which the mobile device and 
base station communicates. The selected channel is used exclusively by the 
specific mobile device and base station until the communication link is terminated 
by either device. No other device can use the channel while the communication 
link is established. 

[0004] In some SISO systems, the selection of the channel is based on the 
strength of the signal received by the base station. By selecting the channel with 
the strongest signal, the bit rate of the communication is maximized. In one 
scenario, the several base stations may receive a signal from a particular mobile 
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device and the base station with the strongest signal will establish a channel with 
the mobile device. 

[0005] Multiple-Input, Multiple-Output (MIMO) systems, in contrast, are 
characterized by base stations and mobile devices having more than one 
antenna. In a typical situation, the transmitted data stream is demultiplexed into 
n T separate sub-streams. Each sub-stream is encoded into channel symbols 
and transmitted over one of the n T transmission antennas. The advantage of the 
MIMO system over the SISO system is the linear increase in spectral efficiency 
for MIMO systems compared to the logarithmic increase in spectral efficiency for 
SISO systems as n T increases. 

[0006] The single-user MIMO link may be described by the equation: 

y = Hx + n (1) 

where y is a n R x 1 receive vector, x is a n T x 1 transmit vector, n is a n R x 1 
additive white Gaussian noise (AWGN) vector, and H is a n R x n T channel matrix. 
An element in the channel matrix is denoted by hy and represents the complex 
gain of the channel between the j-th transmitter and the l-th receiver. In most 
situations, H may be assumed to be random and memoryless. It is also 
assumed that the receiver has perfect channel knowledge. The maximum 
channel capacity is attained when the transmitter also has perfect channel 
knowledge and transmits over k antennas, where k is the rank of the channel 
matrix. 

[0007] The requirement of perfect channel knowledge cannot be met with 
current deployable systems. Therefore, there is a need for MIMO systems that 
can operate efficiently under a feedback capacity constraint. 

SUMMARY 

[0008] One embodiment of the present invention is directed to a MIMO 
system comprising: at least one mobile unit configured to transmit a signal based 
on a selected power matrix and a selected beam-forming matrix; and a base 

-2- 



station configured to receive the signal from the at least one mobile unit, the 
base station comprising: a CSIR estimator configured to receive the signal and 
estimate a channel matrix for the at least one mobile device; and a CSIT 
generator configured to broadcast an index to the at least one mobile unit, the 
index corresponding to a region of a channel matrix space enclosing the 
estimated channel matrix of the at least one mobile unit; wherein the selected 
power matrix and the selected beam-forming matrix are selected based on the 
index broadcast by the CSIT generator. 

[0009] Another embodiment of the present invention is directed to a method 
of optimizing the capacity of a MIMO system, the MIMO system characterized by 
K mobile units, each mobile unit characterized by Q levels, the K mobile units 
transmitting to a base station having n T antennas configured to receive a signal 
transmitted by the K mobile units, the method comprising: (a) generating at least 
one positive, semidefinite channel matrices of size n T x n T ; (b) initializing at least 
one power control matrix and at least one beam-forming matrix; (c) partitioning a 
K x n T x n T channel matrix space into a set of partition regions; (d) generating at 
least one power control matrix and at least one beam-forming matrix based on 
the at least one channel matrix, the at least one power control matrix, the at least 
one beam-forming matrix, and the set of partition regions; (e) repeating steps c 
and d a preselected number, N g , of times; (f) estimating a capacity of the at least 
one generated power control matrix, the at least one generated beam-forming 
matrix, and the set of partition regions; (g) repeating steps a through f a 
preselected number, N s , of times; and (h) selecting the at least one generated 
power control matrix, the at least one generated beam-forming matrix, and the 
set of partition regions corresponding to the largest capacity. 

[0010] Another embodiment of the present invention is directed to a MIMO 
wireless base station, comprising: at least one antenna configured to receive a 
signal from at least one mobile unit; a receiver configured to receive the signal 
from the at least one antenna; a CSIR estimator configured to estimate a channel 
matrix for the at least one mobile unit based in part on the signal from the at least 
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one antenna; and a CSIT generator configured to broadcast to the at least one 
mobile unit a partition index corresponding to a partition containing the channel 
matrix of the at least one mobile unit. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0011] The invention will be described by reference to the preferred and 
alternative embodiments thereof in conjunction with the drawings in which: 

[0012] Fig. 1 is a block schematic diagram of one embodiment of the present 
invention; and 

[001 3] Fig. 2 is a flow diagram of another embodiment of the present 
invention. 

DETAILED DESCRIPTION OF THE PREFERRED AND ALTERNATIVE 

EMBODIMENTS 

[0014] As used hereinafter, a MIMO multi-access system includes K users 
and a single base station in a wireless network. A single base station is used to 
simplify the description of the invention and generalization to systems including 
more than one base station is within the scope of the present invention. Fig. 1 is 
a block schematic diagram of one embodiment of the present invention. For 
purposes of clarity, only two mobile units (users) 101 102, labeled 1 and K, 
respectively, are shown in Fig. 1 , but it should be understood that there are K 
users in communication with the base station. As used hereinafter, user k refers 
to any one of the K users in the MIMO system. 

[0015] Referring to Fig. 1 , the base station 1 10 is equipped with n R antennas 
112 and are referred hereinafter as the receive antennas although it should be 
understood that antennas 112 both transmit and receive signals to and from the 
users. The receive antennas 1 12 are termed receive antennas using a 
convention wherein the base station receives a signal transmitted by the user. 
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[0016] The receive antennas 112 receive the signals, Y, transmitted by the K 
users, indicated in Fig. 1 by 165 for user 1 and 167 for user K, and pass the 
signals 165 167 to the receiver 1 14 and the Channel State Information - 
Receiver (CSIR) estimator 116. The receiver 1 14 separates the received signal 
into the K user data streams based on the channel information provided by the 
CSIR estimator 116 and sends each user data stream to a decoder 118 where 
the channel symbols are decoded to reconstitute the user's data stream. The 
implementation and operation of the receiver 1 14 and decoder 1 18 are known to 
one of skill in the art and examples of such implementations may be found in 
Tran, T.A.; Sesay, A.B, "A generalized simplified ML decoder of orthogonal 
space-time block code for wireless communications over time-selective fading 
channels," Vehicular Technology Conference, 2002. Proceedings. VTC 2002- 
Fall, 2002 IEEE 56th , Volume: 3 , 2002 Page(s): 1911-1915 vol.3 and is herein 
incorporated by reference. 

[0017] The CSIR estimator 116 estimates the channel matrices for each user 
and provides the channel matrices, {vi, . . Vk} to the receiver 114 and Channel 
State Information - Transmitter (CSIT) generator 1 20. In one embodiment of the 
present invention, the CSIR for user k may be obtained by the following equation. 

v,=[y,(l)....y,(M^ (2) 

where x k (m) is the n T x1 m-th transmitted pilot preamble vector of user k, yk(m) 
is the n R x 1 m-th received pilot symbols where M > n T . Alternative methods for 
estimating v k known to one of skill in the art are intended to be encompassed 
within the scope of the invention. For example, Barroso, V., Xavier, J., "Blind 
identification ofMIMO channels: a closed form solution based on second order 
statistics" Signals, Systems, and Computers, 1999. Conference Record of the 
Thirty-Third Asilomar Conference on , Volume: 1 , pp 70 - 74, 1999 discloses 
alternative methods for channel state estimation and is herein incorporated by 
reference. 
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[0018] Equation 2 shows that v k is an n R x n T dimension channel matrix 
characterizing the channel between the k-th user and the base station 110. In 
some embodiments of the present invention, the channel is modeled as a block- 
fading channel where the encoding frame and decoding frame spans across 
many fading blocks. The channel matrix, v k , is highly correlated within a fading 
block any may become less correlated between fading blocks. The channel 
model may be applied to model a slow mobility frequency-hopping channel. 

[0019] The CSIT (Channel State Information - Transmitter) generator 120 
partitions the K x n R x n T dimension channel matrix into Q K regions where Q =2 C 
is the number of quantization levels per user allowed in the system and c is the 
maximum number of bits allowed to feedback per user. The number of 
quantization levels per user represents the number of available power control 
and beam-forming matrices that each user may select from in response to the 
feedback signal received from the base station. Each partition, R q , is assigned a 
label, q, having a value from 0 to Q K -1 . 

[0020] The CSIT generator 120 determines the region containing the v k 
estimated by the CSIR estimator 116 and assigns the corresponding partition 
index, q, to the user's feedback signal, U k . In a preferred embodiment, the 
feedback signal, U k = q k where q k is the value in the k-th position of the partition 
index. The set of user feedback signals 125, (Ui, U 2l . - ., U K ), is broadcast by 
the base unit 1 10 to all K users. 

[0021] On the user side, the index decoder 152 receives the feedback signal 
125 from the base unit 110 and extracts the user's partition index from the 
feedback signal 125. The user's partition index, indicated by 153 for user 1 and 
by 157 for user K in Fig. 1 , is provided to the user's power & beam-forming 
control 154. The user's power & beam-forming control 154 may generate and 
send the power control matrix, pi(qj), and the beam-forming matrix, B r (qj), both 
indicated by 155 for user 1 and by 159 for user K in Fig. 1, to the transmitter 158. 
In some embodiments, the set of pj(qi) and Bj(q { ) are pre-computed and stored in 
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the power & beam-forming control 154. The power & beam-forming control 154 
retrieves and sends to the transmitter 1 58 the pair of pi(q,) and Bi(qj) 
corresponding to the received partition index for the user. 

[0022] In some preferred embodiments, pi(qO and Bj(qO are pre-computed by 
the base station based on the number of users, K, in the base station's cell. It 
should be understood that not all K users transmit their payload simultaneously 
because some of the K users will be in an idle or waiting state that does not 
require transmission. Up to N users may transmit their payload simultaneously 
and in many situations, it is expected that K >N. If another mobile unit enters 
the base station's cell, the base station computes a new set of pi(q,) and Bj(qj) for 
K+1 users and transmits the new pj(qi) and Bi(q,) to each of the K+1 mobile units. 
The storage requirement for the transmission sets, pj(q,) and Bj(qi), remains 
independent of K. 

[0023] The user's data stream is encoded by the encoder 156 and the 
transmitter 158 transmits the encoded data stream over the n R antennas 160 
based on the user's pjfa) and Bj(qj). The transmitted signal, indicated by 165 for 
user 1 and by 167 for user K in Fig. 1 is transmitted to the base station 110. 

[0024] The transmission scheme of all the K mobile users are completely 
specified by {(v^o,...,^^)):^ e[o,e-i]} , the set of Q K partitions of the channel 

matrix space generated at the base station, and by {( Pl ),...., p K (g K )):q k e[o,e~i]} 

and {(b, (q x ),...., (q K )) : q k e[O,0-i]} , which are the Q K sets of power control 

matrices and beam forming matrices, respectively, stored at each of the K mobile 
transmitters. 

[0025] A method for generating the optimum set of parameters is now 
described with reference to Fig. 2, which is a flow diagram of one embodiment of 
the present invention. The parameters that are optimized are the Q K partitions, 
R q , of the K x n R x n T channel matrix space denoted as 
{(v, (<?,)>• .*,v K (q K )):q k e [0,2-1]} , the Q K sets of power control matrices of the K 
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mobile users denoted as {(a(^i)--Pk(^))-^ e[0,e-i]} , and the Q K sets of beam- 
forming matrices denoted as {{B 1 (q l ) 9 .... 9 B K (q K )):q k e[o,e-i]}. Although the 

partitions are stored at the base station and the power control and beam-forming 
matrices are stored at each of the K mobile users, the partitions, power control 
matrices, and beam-forming matrices are generated together iteratively by the 
method. 

[0026] A partition index vector, denoted by q, may be expressed as q = 
(qi,...,qi<) where the k-th component q k e[0,Q-1] is the k-th digit of q expressed in 
base-Q. As an illustrative example, if number of quantization levels per user is 2 
and the number of users is 4, then Q = 2 and K = 4 and a partition index of 5 may 
be expressed as q = 5 = (0,1,0,1) where qi=1, q 2 =0, q 3 =1, and q 4 =0. 

[0027] A distance measure is defined by equation (2) below: 
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where (v-i v K ) denotes the K channel matrices wherein each channel matrix is 

an n R x n T dimension matrix, I k „t is the identity matrix of size Kn T , A* is a system 
parameter related to the average power constraint of user k, and a z is the 
channel noise power. It should be understood that the present invention is not 
limited by the distance measure described in equation (2) and any suitable 
measure apparent to one of skill in the art may be used and is within the scope of 
the present invention. 

[0028] The distance measure is used to partition the K x n R x n T space. The 
q-th partition region, Rq, is defined as the region of channel matrices nearest, in 
the sense of the distance measure, to (p ?i ,..,p 9/t ;B ?i ,...,B 9jc )and is given by 
equation (3) below: 
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[0029] Referring to Fig. 2, a set of Q T partition regions are initialized in step 
210 where Q T = Q K . In step 210, a set of Q positive semidefinite matrices, 
denoted by {A(qj), qje[0,Q-1]} are randomly generated in step 210. Each matrix, 
A(qj), is of size n T x n T . The set of power control matrices for user k in a partition 
qi is denoted by p k (qO and consists of the matrices {p k (0),..., p k (Q-1)} wherein 
each pk(qi) corresponds to the eigenvalue matrix of A(q k ). The power control 
matrices are diagonal matrices of size n T x n T having diagonal elements equal to 
the n T eigenvalues of each A(q k ). 

[0030] The set of beam-forming matrices for partition q\ for user k is denoted 
by B k (qj) and comprises the set, {B k (0),...,B k (Q-1)}. Each B k (qj) is initialized as 
the hermitian matrix corresponding to the eigenmatrix of A(qi) for all qi=[0 f Q-1]. 
The eigenvalues and eigenmatrix of A(qj) may be obtained using methods known 
to one of skill in the art such as, for example, single value decomposition and do 
not require further discussion. The partitions are generated using equation (3) 
and the Q matrices, A(qj). 

[0031] Once the partitions, {Ro,...,Rqt-i}, are initialized in 210, in step 220, a 
new set of power control and beam-forming matrices are calculated by iteratively 
solving the system of equations below: 



(p 9 ,B 9 ) = argmin 



Z ^(v 1 ,.,v j ,)^(^(( v i 5 - v k);(a ? -- 5 p^);(Bi,-.,b /c ))) 

qeO(q k ) 



(5) 



[0032] The new set of power control and beam-forming matrices obtained in 
step 220 are used to define a new set of partitions in step 230 by applying 
equation (4) to the new set of power control and beam-forming matrices. 
Convergence of the matrices and partitions are tested in step 235. If the 
matrices and partitions have not converged, steps 220 and 230 are repeated 
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before the convergence is tested again in step 235. In a preferred embodiment, 
steps 220 and 230 are repeated N g times where N g is a pre-selected value. 

[0033] Step 240 estimates the capacity of the system with the converged 
power control and beam-forming matrices from step 220 and the partition regions 
from step 230. In a preferred embodiment, the capacity is estimated by the 
equation below: 

q=Q 

[0034] Steps 210 through 240 may be repeated for different initializations of 
{A(qj)}. Each initialization produces a set of power control matrices, beam- 
forming matrices, and partition regions along with an associated capacity. Step 
245 branches to step 210 if a new initialization is desired or branches to step 250 
if a sufficient number of initializations have been performed. The number of 
initializations may be a pre-selected value and may be determined by one of skill 
in the art by considering and balancing factors such as, for example, desired 
sample coverage and computational cost. 

[0035] The set of power control matrices, beam-forming matrices and partition 
regions having the largest capacity is saved in step 250. In some preferred 
embodiments, the power control matrices and beam-forming matrices are 
downloaded to each mobile device. 

[0036] Having described at least illustrative embodiments of the invention, 
various modifications and improvements will readily occur to those skilled in the 
art and are intended to be within the scope of the invention. Accordingly, the 
foregoing description is by way of example only and is not intended as limiting. 
The invention is limited only as defined in the following claims and the 
equivalents thereto. 



log 2 



0 . 0 
0 0 p K (q K )) 



(5) 



- 10- 



